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FORTRAN PROGRAM FOR COMPUTATION OF GROUP 


TABLES - ALPHANUMERIC DISPLAY 
by Gabriel Allen and David D. Evans 
Lewis Research Center 

SUMMARY 

An existing FORTRAN program for machine computation of finite groups has been 
modified to allow the output to be displayed as symbols rather than as numerical cycles. 
The program is written for second generation machines. As written, the program 
makes use of, but does not depend on, the fact that every finite group G of order n is 
isomorphic to some subgroup Tj of the symmetric group S fi . The program can print 
a table of any possible combination of subsets of group elements of G. The procedure 
for using the program is as follows: After finding the Tj, those elements of Tj which 
correspond to the subsets of G whose table is desired are read in as input data in cycle 
notation. The symbol for each element of G as it is to appear in the table is also read 
in. The group multiplication table for the subsets of G is computed and printed using 
these symbols. Examples are shown of complete group tables, tables of the group 
multiplication of subsets and a change in symbol notation. 


INTRODUCTION 

Since its first timid beginnings sixty years ago, group theory has become a standard 
branch of applied mathematics. Reference 1, which was published in 1959, has a list 
of 213 references of applications of group theory. This list (which is not exhaustive) 
includes applications to crystallography, chemistry, atomic and molecular spectro- 
scopy, elementary particles, and relativity. However, the computational effort in- 
volved in dealing with large finite groups can be prohibitive. As a result, computing 
machines have been applied to the task of proving theorems and algorithms about finite 
groups whose group tables were entered into the program as input (ref. 2). 



Recently, a FORTRAN program was written for the machine computation of the 
group table itself, using only the elements of the group as input data (ref. 3). This 
program was also made available for use on second generation machines (ref. 4). The 
method made use of the fact that every finite group of order n is isomorphic to some 
subgroup of the symmetric group S n (ref. 1, p. 260). Therefore, in using this pro- 
gram, the input data were entered in the form of cycles, the elements of S n » The out- 
put was also printed in the form of cycles and the isomorphism between the elements of 
G and S n was then used to construct the desired table for G. This last step had to be 
performed manually. 

This program saved the bulk of the effort required to construct a group table. 
Nevertheless, the fact that the entries in the printed output appeared in the form of 
cycles meant that a decoding process (the isomorphism between G and the appropriate 
subgroup of S n ) had to be used to write the final table for G. This decoding process 
could still be time consuming. The modified program (which is written for second 
generation machines) performs this decoding process and prints the group table in 
easily recognized alphanumeric symbols rather than numerical cycles. 

The major new function of the program is contained in a subroutine with two entry 
points. This subroutine identifies a given element of Tj and substitutes for it the 
symbol for the appropriate element of G in the output format. 

With the use of this program, a desired group multiplication table can be obtained 
in the following manner: 

(1) A set of convenient alphanumeric symbols for the elements of G is selected. 

(2) An isomorphism between G and some subgroup T. of S. is found. Here S. 

3 J 3 

is the smallest symmetric group which can serve this purpose. 

(3) The elements of Tj which correspond to the elements in the subsets of G 
whose multiplication table is desired are entered as input in the form of cycles. 

(4) The symbols chosen in step 1 are entered as input. 

(5) The group multiplication table is computed and the output printed using the 
symbols chosen in step 1. 

An example of the group multiplication table of two subsets of a group is presented 
in the output. 


SYMBOLS 

A n alternating subgroup of symmetric group, S n 

C k j element in row and 1 th column of matrix C 

Dg dihedral group 

G finite group of order n 
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O the cubic group 

S n symmetric group of order n- factorial 

Tj some subgroup of S n which is isomorphic to O 
twofold rotation axes 


DEFINITIONS, CONVENTIONS AND PREPROGRAM PREPARATIONS 

The symbol S n designates the permutation group of n objects and is of order n 
factorial. The even permutation group of n objects A n is a proper subgroup of S n . 
The numbers between commas in a cycle will be called units. The convention adopted 
herein for describing the effect of a cycle is that in which the units denote objects and 
each unit is moved to the location currently occupied by the unit to its left. 

Before using the program, some suitable alphanumeric symbols for the elements 
of G are selected. These symbols would preferably be as close to the symbols in the 
literature for that group as the printer is capable of printing. Next, an isomorphism is 
found between G and T.., the subgroup of some symmetric group Sj, which is as small 
as can be thought of for the purpose. Finally, the particular sequence of units in a 
cycle for a given element of must be decided upon, for the main program must store 
the cycles in this exact form in the doubly dimensioned array representing the group 
table. See appendix A for an example of how this can be done for S^. 

The user normally has the freedom to choose any sequence of units for an element 
of Sj which forms a cycle equivalent to the one used herein (see ref. 3). However, in 
this program the smallest integer in the entire cycle must always be the left- most one. 

It should also be noted that the largest group which this program will handle is S^. 
Once the principles of the particular decoding process used in this report are under- 
stood by the user, the extension to larger groups should be straightforward. However, 
such an extension will require that subroutine READ be rewritten. 


PROGRAM DESCRIPTION 
General Description 

After the group multiplication has been performed, the results are stored in the 
doubly dimensioned array C. An entry in this array C M corresponds to the element 
in the k in row and the l in column of the table. However, in this array, the element 
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is stored as a cycle chain. The program identifies Cj^ as a specific member of 
and prints in its place that alphanumeric symbol chosen to represent the element of G 
corresponding to the identified element of S^. This process will be clarified subse- 
quently by an example. 

The program as presented in this report makes use of the symmetric group to ob- 
tain a group table for Tj which remains internal to the machine. The decoding sub- 
routine can be applied to the translation of any list of cycles as input to any list of de- 
sired symbols. All that is necessary is a correspondence between the desired symbols 
and the cycles . As a matter of fact, some of the complications of the program in ref- 
erence 4 were needed merely to simplify the output format. In this report, those parts 
of the program have been eliminated, resulting in a simplification. 

Those parts of the program needed for understanding its use are described next. 
For more detailed information, reference may be made to the listing in appendix B and 
to the flow chart of subroutine READ in figure 1. The examples are contained in a 
separate section following this main section. 


Main Program 

Block 1 - Set up constants . - The constants used in the program are given literal 
names and are declared either INTEGER or LOGICAL in TYPE statements. All the 
subscripted variables are dimensioned and alloted storage locations by the use of 
DIMENSION statements. The names of the variables of the group elements of are 
equivalenced to members of the dimensioned variable INPT. A DATA statement is used 
to give literal names to the following variables : 


Variable 

Literal name 

Identification 

( 

LP 

Left parenthesis 

) 

RP 

Right parenthesis 

J 

CM 

Comma 

, 

PER 

Period 


BLK 

Blank 

E 

IDNT 

Identity element of group 


Block 2 - Read in data . - The data describing the group, the group elements, and 
the decoding process are read into the program. The group itself is defined by the 
variable XDENT. The assignment of labels to the objects and to the locations are read 
in from the list for IDENT which is limited to one card. No shifting or ’SQUEZ'ing is 
required, thus simplifying the programming compared to reference 4. The input cards 
are read in the following order: 
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First card: TITLE, one card with FORMAT (80A1) 
Second card: IDENT, one card with FORMAT (80A1) 


The input for IDENT must be long enough to include the maximum number of objects to 

be moved in any group operation. Thus, the list for IDENT must be of the form(Ap 
A 0 , . . . , A.) where j is the subscript of S, and A, is any alphanumeric symbol. 

Several error checks are run on IDENT. These include checking the maximum 
number of nonblank units (which is six in this version), the presence of blank spaces on 
the card, and the illegal use of parentheses. 

The set of group elements whose multiplication table is desired (this may be a sub- 
set of the entire group) is read in the next two sets of cards. These two sets are 
identified in the program as INPUT(J,K). The first of these sets contains the ACROSS 
(column) operations and consists of the INPUT cards for which K = 1. There can be as 
many as four cards with FORMAT(80A1), and each operation is followed by a period. 

A blank card follows the preceding set. It is needed to signal the end of the ACROSS input. 

The next set contains the DOWN (row) operations and consists of INPUT cards for 
which K = 2. There can be as many as four cards with FORMAT(80A1), and again, 
each operation is followed by a period. 

A blank card follows the preceding set, signaling the end of the DOWN input. 

The input for both ACROSS and DOWN operations is in cycle notation with the 
smallest integer consistently to the left. No blanks are allowed between the first left 
parenthesis and the final right parenthesis in any cycle chain. If integers are not used 
as the units in the cycle, subroutine READ will have to be rewritten to identify the units 
which are used. The complete set of column operations is read first. Then the com- 
plete set of row operations is read. Each of these sets is limited to four cards and is 
followed by a blank card. 

The FORMAT for the input list is A6. All entries must be left- adjusted in order 
that the output be in line. The example of input which follows is a class product block 
from (see ref. 5). It is the product of the class Rl, R2, R3 (see appendix A) times 
the class containing the rotations about all the threefold axes. 


Card 


1 |S4 || 

2 

3 |' % £ ? 4 ? J3p . K i ? £ « 3 ) , l > 1 ? 4 ? £ ) . k l ? 3? 4 ) 7 ( 2V3 ? 4 ) . 7 1 ? 3 ? £ ) • (, 1 * £ j 4} » (. 1 ? 473 ) I 

4 


5 

|UV3?ppv4}:U5£> 

{ o ? 4 ) 7 i 5 4 y £ s 3 ) * 
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7 

|I 1 | |I2 13 

14 15 16 

SRI 

SR£ 

SR3 

SR4 

SKI 

8 

ISR3S | ER4S Rl 

R2 R3 RiS 

R2S 

R3S 

RiC 

R£L 

“R3C 


I 
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Card 3 contains the ACROSS (column) operations, card 5 contains the DOWN (row) 
operations, and cards 7 and 8 contain the list of ' ’values” for ENPT. Note that the en- 
tire group is listed even though only two classes are multiplied. This is done because 
classes are not closed with respect to group multiplication and elements of which 
are not in either of the classes used as input can appear in the product. 

See reference 4 for a description of the rest of the main program. 

The data for the decoding process are read from the list for the variable INPT. 

The list for INPT consists of the actual names desired for the group elements in the 
output. The list, of course, assigns ” values” to the dimensioned variable INPT. The 
’’variable” names are listed in an equivalence statement in block 1. In this program, 
the alphanumeric data in the list for INPT (the values ) are the same as the symbols in 
the equivalence statement (the variables). For example, the 11 th value on the list for 
INPT is SR1S; and in block 1, INPT(ll) is equivalenced to the variable SR1S. 

This procedure makes it convenient to use any other choice of variables which are 
isomorphic to the group in place of the ones used here. If the symbol Q5 were de- 
sired in place of SR1S, assigning the value Q5 to the variable SR1S would cause Q5 
to be printed everywhere that SR1S appears now. The important point is that no change 
in the body of the program itself is required. The single change in the input list takes 
care of the entire new decoding process. 

An example of such a substitution is shown in the output for Sg. In the example, 
the well-known isomorphism between Sg and the dihedral group Dg is used to provide 
symbols for the elements of Sg which are more covenient than the ones which arise 
when Sg is regarded as a subgroup of S^. The common symbols for the elements of 
Dg are used. These are a threefold rotation axis R and three twofold rotation axes 
Pp Pg, and Pg. This new output was obtained by simply changing cards 7 and 8 to read 
as follows: 

IRH01 |I£ 13 RH03 15 RHD2 SRI R SR3 SR4 "SR IS RS j 

SR3S |SR4S R1 R£ R3 R1S R2S R3S R1C R2C R3C | 

Subroutine READ 

This is the part of the program which decodes the cycles stored in INPUT and OUT 
(see listing) into the desired symbols for the elements of G. A process of elimination 
is used to identify the cycle, character by character. Advantage is taken of the fact 
that all equivalent cycles are stored in the same unique form (see ref. 3). Each integer 
in a given cycle is identified in succession beginning with the left- most one. The identi- 
fication of an integer is a simple comparison process between the values assigned to the 
variables ’’ONE, ” "TWO, ” ’’THREE" and the unit of the cycle. When needed, left 
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and right parentheses are identified by comparing them against the values of the vari- 
ables LP and RP, respectively. After a given character has been identified, control 
is shifted to another set of comparison tests to be made on the next character in the 
cycle. As each character of the cycle is identified, the number of possible group ele- 
ments to which the cycle could correspond decreases. The identification is complete 
when there is only one group element which could correspond to the given cycle. 

To keep the program as short as possible, it is assumed that certain errors will 
not be made. For example, if 4 is the largest integer to appear in a cycle, tests are 
made only up to and including the integer 3. If the number in the cycle fails all of the 
tests , it is assumed that it is a 4 in the identifying process. 

Illustrative example 1 - Identification of the cycle (1, 3, 4, 2) as the group element 
Rg„ - The cycle (1, 3, 4, 2) enters the subroutine as the value of the variable INPUT in 

the calling vector. The identification proceeds by examining individual characters in 
the array in sequence. This array has four integers , three commas, a left parenthesis, 
and a right parenthesis, for a total of nine characters. The first character to be ex- 
amined is the left- most one. Since it is a left parenthesis, it fails all the comparison 
tests and the next character is examined. This one passes the comparison test on the 
integer 1 and control is thereby shifted to statement 1100. The comparison test in this 
series of statements takes into account the fact that some unit has been identified before 
transfer to statement 1100. Since the first integer in any cycle is always followed by a 
comma, the subroutine skips the next entry and goes to entry 1 + 2 for the next com- 
parison test. Since this entry is the integer 3, control is shifted to statement 1130. 

The (I + 3) rc * entry in the array is examined next. It is not a right parenthesis , 
so control is shifted to statement 1140. The (I + 4) th and (I + 5) th entries both fail 
their comparison tests and the array is thus uniquely identified as the group element Rg< 

Illustrative example 2 - Identification of the cycle chain (1, 3)(2, 4) as the element 
R1S . - This particular array has a total of 10 characters divided amongst four inte- 
gers, and two each of commas, and left and right parentheses. The first entry in the 
th 

array (the K one) is again a left parenthesis . After it has been identified, the 
(K + 2) n< ^ entry (K + 2) is examined. Since it is a 1, control is shifted to statement 
1100. As in the previous example, the next entry to be examined is the (K + 2) n< ^ one. 
Since it is a 3, control is shifted to statement 1130. After the comparison tests in this 
series are completed, control is shifted to statement 1131, where the identification R1S 
is made. 
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EXAMPLES OF OUTPUT 


Full Group 


S4 

« 1,2,3, A) 

( H2). (3, 4 ).( 1,41.12,3). (2*41. <1,3 >.(2,4,3 ).( 1,2,3) . ( 1,4,2 ).(1,3.4K 
(2*3, 4 1. <1*3,21.(11 2, 4). <H4, 3). <1,2, 3, 4). (1,4,2,31. 11, 3, 4, 2). (1,3)12,4). 
<1121 (344) .(1,4) <2131.(1,4*3,23. <1,3,2, 4). <1,2, 4,3). 


( 1*2 ). J 3, 4 }.(1„ 4). 12,3). (2*4) .11,3)- (2.4 « 3 ) . ( 1,2,3) . < 1,4,2) . (1,3,4) . 
<213,4).<U3,2).U42,4).( 1,4,3). <1,2,3,4). ( l,4,2,3).< 1,3,4,2>. < 1,3)(2,4). 
(E12) (3,4 I • ( 1 » 4 M2 »3>. <1* 4b 3*29. ( 1, 3,2,4) . U *2,4,3) - 


S4 



11 

12 

13 

14 

15 

16 

SRI 

SR2 

SR3 

SR4 

SR1S 

SR2S 

SR3S 

SR4S 

RI 

R2 

R3 

RIS 

R2S 

R3S 

11 

E 

R2S 

SR3 

SR2 

SR3S 

SR2S 

R3C 

14 

13 

R3 

RI 

16 

15 

RIC 

SR1S 

R3S 

SR4 

R2C 

12 

R2 

12 

R2S 

E 

SR4 

SRI 

SR IS 

SR4S 

14 

R3C 

R3 

13 

15 

RIC 

Ri 

16 

SR3S 

RIS 

SR3 

R2 

11 

R2C 

13 

SR 3 S 

SR4S 

E 

R3S 

SR3 

SR4 

RIC 

RI 

!*? 

16 

R2 

R2C 

11 

12 

SR2 

SR IS 

RIS 

R 3 

R3C 

14 

14 

SR2S 

SR1S 

R3S 

E 

SRI 

SR2 

15 

16 

RIC 

RI 

12 

11 

R2C 

R2 

SR4 

SR4S 

R2S 

R3C 

R3 

13 

15 

SR3 

SRI 

SR3S 

SR 1 S 

E 

RIS 

12 

R 2 

11 

R2C 

14 

R3 

13 

R3C 

R3S 

SR2 

SR2S 

16 

ri: 

RI 

16 

SR2 

SR4 

SR4S 

SR2S 

RiS 

E 

R2C 

11 

R2 

12 

R3 

14 

R3C 

13 

R2S 

SR3 

SRI S 

15 

Ri 

RIC 

SRI 

RIC 

15 

R2C 

12 

14 

R3C 

SR1S 

SR4S 

SR2S 

SR3S 

E 

R2S 

R3S 

RIS 

13 

16 

U 

SR.2 

SR3 

SR4 

SR2 

16 

R 1 

R2 

I 1 

R3C 

14 

SR3S 

SR2S 

S&4S 

SR IS 

R2S 

E 

RIS 

R3S 

R3 

RIC 

12 

SRI 

SR4 

SR3 

6R3 

15 

RIC 

11 

R2 

13 

R3 

SR4S 

SR1S 

SR3S 

SR2S 

R3S 

RIS 

E 

R2S 

14 

RI 

R2C 

SR4 

SRI 

SR2 

SR4 

R 1 

16 

12 

R2C 

R3 

13 

SR2S 

SR3S 

SR IS 

SR4S 

RiS 

R3S 

R2S 

E 

R3C 

15 

R2 

SR3 

SR2 

su 

SR1S 

R 3 

14 

R1 

15 

12 

R2 

E 

RIS 

R2S 

R3S 

SRi 

SR3 

SR4 

SR2 

R2C 

R3C 

RIC 

SR4S 

SR2S 

SR3S 

SR2S 

14 

R3 

RIC 

16 

R2C 

II 

RIS 

E 

R3S 

R2S 

SR4 

SR2 

SRI 

SR3 

12 

13 

RI 

SR3S 

S RI S 

SR4S 

5R3S 

13 

R3C 

15 

R 1 

ll 

R2C 

R2S 

R3S 

6 

RIS 

SR2 

SR4 

SR3 

SRI 

R2 

14 

16 

SR2S 

SR4S 

SR 1 S 

SR4S 

R3C 

13 

16 

RIC 

R2 

12 

R3S 

R2S 

RIS 

E 

SR3 

SRI 

SR2 

SR4 

11 

R3 

15 

SR 1 S 

SR3S 

SR2S 

R 1 

SR4 

SR2 

SRI S 

SR3S 

R2S 

R3S 

11 

R2C 

12 

R2 

R3C 

13 

R3 

14 

RIS 

SRi 

SR4S 

RIC 

16 

15 

R2 

R IS 

R3S 

SR2 

SR3 

SR4S 

SR1S 

13 

R 3 

R3C 

14 

RIC 

15 

16 

RI 

SR2S 

R2S 

SRI 

12 

R2C 

I 1 

R3 

6R1S 

SR2S 

R2S 

R1S 

SR4 

SR3 

16 

15 

RI 

RiC 

R2C 

R2 

12 

11 

SRI 

SR3S 

R3S 

13 

14 

R3C 

R1S 

R2 

R2C 

R3C 

R3 

16 

15 

SR4 

SR3 

SR 2 

SRI 

SR2S 

SR IS 

SR4S 

SR3S 

RIC 

11 

14 

E 

R3S 

RZS 

R2S 

12 

I 1 

R3 

R3C 

Ri 

RIC 

SR2 

SRI 

SR4 

SR3 

SR3S 

SR4S 

SR1S 

SR2S 

15 

R2C 

13 

R3S 

E 

RIS 

R3S 

R2C 

R2 

14 

13 

RIC 

Ri 

SR3 

SR4 

SRI 

SR2 

SR4S 

SR3S 

SR2S 

SR IS 

16 

12 

R3C 

R2S 

RIS 

E 

RIC 

SRi 

SR9 

SR2S 

SR4S 

R3S 

R2S 

R2 

12 

R2C 

11 

13 

R3C 

14 

R3 

E 

SR4 

SR3S 

RI 

15 

16 

R2C 

R3S 

R 1 $ 

SRI 

SR4 

SR2S 

SR3S 

R3 

13 

14 

R3C 

16 

Ri 

RIC 

15 

SR4S 

E 

SR2 

I 1 

R2 

12 

R3G 

SR4S 

SR3S 

R IS 

R2S 

SR2 

SRI 

RI 

RIC 

16 

15 

II 

12 

R2 

R2C 

SR3 

SR2S 

E 

14 

13 

R3 



RIC 

R2C 

II 

SR4S 

RIS 

12 

SR2S 

R3S 

13 

SRI 

SR2S 

14 

SR3 

SRBS 

15 

R2S 

SR4 

16 

R3S 

SRI 

S«1 

R 3 

RI 

SR 2 

13 

15 

SR3 

R3C 

16 

SR4 

14 

RIC 

SRI’S 

l l 

13 

5R2S 

R2 

R3C 

SR 3 6 

12 

R3 

SR4S 

R2C 

14 

RI 

e 

SR3 

R2 

SR3S 

E 

R3 

SR2 

SR4S 

RIS 

RI 

12 

R2S 

16 

R2 

R3S 

15 

1 1 

RIC 

RIS 

SR2 

R2C 

SRI S 

R2S 

R3C 

SR4 

SR 1 S 


R3C 


SRI 
SR2 
R2S 
R1S 
SR4S 
SR3S 
R2 
R2C 
12 
! 1 
16 
15 
RiC 
R1 

SR2S 

SR4 

E 

r3 

14 

R3 

SR1$ 

SR3 

R3S 
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Class Algebra Block of 


S5 

i H2,‘3»5> 

<2*4»3).(H2,3a.(W4,2>.(H3,55. (2,3 t 4MI,3,2I.J l * 2,4 ) . (1 ,4, 3 ) . 
(li3U2*4).U,2M3*4)« <1,4H2,3)« 


S4 



SRI 

SR2 

. « - . X .. 

SR3 

SR4 

SR IS 

SR2S 

SR3S 

SR4S 

R1S 

R2S 

R36 

SR4 

SR2 

SR3 

SR3 

SRI 

SR4 

SR2 

SR4 

SRI 

SRI 

SR3 

SR2 

SR2S 

SR3S 

SR4S 

SR1S 

SR4S 

SR3S 

SR4S 
SR IS 
SR2S 

SR3S 

SR2S 

SRIS 


Sg as a Subgroup of S^ 


S3 

( 1 * 2,31 

<n2)..ll,3)«t2t3>.4l,2,3). tl f 3t2). 
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SUMMARY OF PROCEDURE 


A FORTRAN program for machine computation of finite groups has been written 
which allows the output to be displayed in alphanumeric symbols rather than in numer- 
ical cycles. The program is written for second generation machines. As written, it 

makes use of the fact that every finite group G of order n is isomorphic to some sub- 
group Tj of the symmetric group S n . However, the method is not dependent on this 
fact. The program can print a group multiplication table of the product of any two sub- 
sets of G. 

The main new feature of the program is an efficient decoding process from cycles 
to any chosen alphanumeric symbols. This decoding process is not limited to elements 
of groups . It can be applied to the translation of any list of cycles as input to any de- 
sired symbols as output. 

The procedure for using the program is as follows : 

1. A set of desired alphanumeric symbols for the elements of G is selected. 

2. An isomorphism between G and some subgroup T n of S n is found. 

3. The elements in the subsets of G whose group multiplication table is desired 
are entered as input in the form of cycles. 

4. The symbols chosen in step 1 are entered as input. 

5. The group multiplication table is computed and the output printed using the 
symbols chosen in step 1. 

The program is complete for groups up to those which are isomorphic to or 
any of its subgroups . A straightforward modification can extend it to larger groups 
if desired. 

A change in choice of symbols for elements of G is very easy (see step 4 in the 
procedure) and an example of such a change is shown. Examples of both a complete 
group table and a group multiplication table of two subsets of a group are also shown. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, October 7, 1970, 

129-03. 
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APPENDIX A 


FORTRAN NOTATION FOR GROUP O FOR CODING 
OF ISOMORPHISM WITH S 4 

In the output examples in this report, the following notation is used: The group O, 
the cubic group, is isomorphic to the symmetric group S 4 . The former group has four 
threefold axes, three fourfold axes and six twofold axes. Since FORTRAN IV does not 
have lower-case letters, the following alphanumeric symbols will be used to represent 
the group elements of O. 

The symbol R denotes a rotation about one of the fourfold axes, SR (for small R) 
represents a rotation about one of the threefold axes, and I represents a rotation about 
a twofold axis. The particular axis in a given case is identified by a number. The 
letter S following the axis identification means squared and signifies that the rotation 
has been repeated about the given axis. The letter C following the axis identification 
means cubed and signifies that the rotation about the given axis has taken place three 
times. Thus, R2C is the group element of O signifying that a rotation has taken place 
three times about fourfold axis number 2; whereas SR1S means that a rotation has oc- 
curred twice about threefold rotation axis number 1. (See ref. 5, p. 27, for a con- 
venient table of group O using this notation. ) 

The following code is used for the isomorphism between cycles (the elements of S^) 
and the elements of O: 

II— (1,2); 12— (3,4); 13 — (1,4); 14— (2,3); 15— (2,4); 16 — (1,3) 

SRI—— (2,4,3); SR2 —(1,2,3); SR3— (1,4,2); SR4 — (1,3,4) 

SR IS — (2,3,4); SR2S — (1, 3,2); SR3S — (1, 2, 4); SR4S —(1,4,3) 

R1S — (1, 3)(2, 4); R2S — (1, 2)(3, 4); R3S— (1, 4)(2, 3) 

Rl — (1,2, 3,4); R2 — (1,4, 2, 3); R3 —(1,3, 4, 2) 

R1C — (1, 4, 3, 2); R2C — (1, 3, 2, 4); R3C— (1, 2, 4, 3) 
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APPENDIX B 


PROGRAM LISTING 


C etOCK(l) SETS UP CONSTANTS 

DIMENSION I TOPI 120) ,IS I D ( 120 ) , I DOT ( 22 ) , L INPUT ( 15), INPTI24) 
DIMENSION INPUT! 400,2) , INITAU6), IDENTI80) , NPAIR(24,2) , 
IL4NE1576),NIN124,2) ,N0PA.(2) ,KN(6) ,N1NPUTI2 ) 

INTEGER RP,CN,TEMP,BfcK, CHAR, PER, PAIRS! 18 , 24 , 2 ) , ANS I 6 , 24, 24 ) , OUT I 15 
1,24 , 24 ) , T I TLE (80 ) , OUTPUT! 21 ) 

INTEGER 6R1,SR2»SR3»SR4»SR1S»SR2S»SR3S»SR4S»R1»R2,R3*R1S»R2S,R3S* 

• R1C »R2C ,R3C 

DATA IDI1/6H1 / 

DATA! IDOT ( l ) J 1 = 1 » 21 )'/21*6H / 

DATA LP,flP,CW,BLK, PER/1 Hi , 1H) , 1H, ,1H , 1H./ 

LOGICAL COMMA, ONE 
£OMMeN/IKYKAX/IKY,MAX,NCHK 
DGMMON/NLMOPS/JO ,K0 ,N0UTI24,24) 

COMMCN/LITRAt/I I 1,112,13414, 15, 1 6, SR 1 , SR2 , SR3 , SR4, SR IS , SR2S, SR3S , 
«£R4S'»R1,-R2,R3,R1S,R2S,R3S,R1C,R2C,R3C 
EQUIVALENCE 4 INP T ( l H 1 1 1 1 

C BLOCK 12) READS IN DATA 

1 READ (5*66) TITLE 
READ (5,66) 1DENT 
MRIT£(6,€0) TITLE 
MR I TE < 6,81) IDENT 
ILP = C 

I RP = C 
N-C 

DHAR*BLK 
00 5 1=1480 
TEMP= IDENT ( I ) 

IF UTEMP.EC.BLK) GO TO 5 
IF (TEMP.EQ.LP) GO TC 2 
'IF (TEMP. EC. RP) GO TO 3 
IF (TEMP. EC. CM) GO TO 4 
OHAR=TEMR 
GO TC 5 

2 ILP=ILP+1 
GO TO 5 

3 IRP=IRP+1 

4 N*=N + 1 

if CN.GT.6i GO TO 61 
IF (CHAR4BC.BLK) GO TO 62 
IN1TAHNHCHAR 
CH AR=BLK 

5 CONTINUE 

IF f ILP.NE. liOR. IRP.NE. 1 ) GO TO 63 
I0-8C 

DO 8 K= 1 4 2 
DO 7 1=145 
US=( l-l )#80+l 
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j£~ I ££0 

READ 45,66) I INPUT C J % K I t J = J $ « JE ) 

MIT£C6 f 90) -4 INPUT I J*K> , J^JStJE) 

00 6 . J* J S f J E 

IF (INPUT! J?K)®NE® BLK I GO TO 7 

6 CONTINUE 

N INPUT I K i^JS 

'GO TO 8 

7 CONTINUE 

8 CONTINUE 
U€=J£^80 

READ! 5, 5678) INPT 
5678 FORMATt 12A6) 

BLOCK 43 ) STORES EACH OPERATION IN PAIRS ARRAY AS A PRODUCT OF 
TRANSPOSITIONS 

DO IS K=a*2 
NPA IP ( 1 *0 
N0P*1 
IPRN-C 
CHAR^BLK 
41*0 

ONE". TRUE. 

COMMA*. FAtSEi 
l€*N INPUT IK) 

DO 16 1*1, IE 
TEMP* INPUT ( I |K) 

IF CTEMP.EQ.8LK) GO TO 16 
IF (TEMP.EQ.fcP) GO TO 9 
IF CTEMP^ EQ ®PP ) GO TO 10 
IF CieMP.EC.CM) GO TO 11 
IF f TEMP. EC. PER) GO TO 14 
CHAR* TEMP 
GO TO 15 

9 IPRN* IPRK+ 1 
ONE"® TRUE ® 

4F ( I PRN*LT • 0*0R . iPRN.GT ® 1 ) GO TO 57 
GO TO 15 

10 4PRN^IPRN-1 

IF ( .NOT wCOMMA) GO TO 58 
COMMA*. FAUSEi 
0NE=i-TRUE* 

IF HPRN.LT.O.OR.IPRN.GT.l) GO TO 57 
GO TO 12 

11 COMMA** TRUE ® 

12 NPA IP (NOP, K )*NPA I R (NOP, KJ+1 

IF (NPA I R ( NO R » K ) « GT * 1 7 ) GO TO 59 
IF fCHAR*EQ s 0LK) GO TO 56 
4 J = NP A I R I NOP | K ) 

PAIR5U J%NOP|K) = CHAR 

IF CONE ) GO TO 13 

NPA IP (NOP ,K)^NPA!R(NCP f Km 

IJ=NPAI RIN0P4K) 

PAIRS II Ji, NOP | K)* CHAR 

13 CHAR*BL K 
0NE^#FALSE® 

GO TO 15 



14 NOP-fvCP+I 

IF (NOP»GT *24} GO TO 60 
NP A I-R ( NOP * K ) +0 

15 II = HU 
INPUT!! I f K } “TEMP 

16 CONTINUE 

N I N P \l 1 i K i n 1 I + l 

If I f NPUT I 1 1 * K) * NE • PER) GO TO 17 

NOP-NOP-1 

GO TO 18 

17 INPUT (I I * 1 *K) =PER 

l l 8 IF IftPAIRU*K).EQ.O) GO TO 1 

NOPAT K ) =NOP 

19 CONTINUE 

C BLOCK (4 ) SETS UP ANS ARRAY IN STANDARD FORM OR CONFIGURATION 

NO P 1=N0P A I 1 ) 

NOP2*NOPAC2) 

DO 2C I 1*1 f NOPl 

00 20 12^1 ? NGP2 
DO 2C I = 1 * N 

20 ANSI !fIHI2)*INITALH) 

C BLOCK i 5 ) PERFORMS THE PRODUCT OPERATIONS 

DO 27 11*1 9 NGP1 
K 1 = NP A I R I 11*1) 

DO 27 12*1* NGP2 
K2=NP AI R TI 2*2) 

DO 2 3 I = 1 ? K 1 1 2 
y* k i™ i+ 1 
y i=o 
y 2 - o 

DO 22 L = 1*N 

IF f ANSI L fill 12) -NE^PAIRSIJ* II* 1 ) ) GO TO 21 

y i = l 

21 IF (ANSlL* 1 1 1 1 2) .NE • PA I RSI y-1 , 1 1 , i ) ) GO TO 22 

y 2 =l 

22 CONTINUE 

IF { ’Jl.EC.O.OR. J2*EQi.O) GO TO 55 

1 T M P = AN S I J 1 * 1 1 * I 2 ) 

IANS { J 1 * II* 12)=AN$<J2%I1* 12) 

ANS C y 2 ? 11* I 2 1 = ITMP 

23 CONTINUE 

DO 26 I S TL * K 2 f 2 

j*k 2 - 1+1 

J1 = 0 
J2 = 0 

DO 25 L = 1 * N 

IF I ANS I L * Il| 12) .NE. PAIRS! J* 12,2)) GO TO 24 

y i=l 

24 IF (ANS(t*Ili 12) .NE. PAIRS! J-1?I2*2) ) GO TO 25 

J2 = L 

25 CONTINUE 

IF fJl*EC®0*GR*J2®EG®0) GO TO 55 
ITMP*ANS! Jl t I It I 2) 

ANSI Jit Ilf I 2 I -AN S C J2 * II * 12) 

■ANSI J2* II* 121 -IT HP 
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26 CONTINUE 

27 CONTINUE 

BLOCK 16 1 IDENTIFIES THE SINGLE GROUP ELEMENT WHICH HAS THE SAME 
EFFECT AS BLOCK! 51 

MA X“C 

DO 38 I l fc I *NGP1 
DO 38 12^1 ? N0P2 
JRP^C 

K 1 = 1 
KA = 1 

OUT ( 1 1 1 L ? I 2 ) ^LP 
DO 28 KK fe 1 ? N 

28 KN i KK ) =KK 
00 35 I ~1 P N 
DO 29 L=1 f N 
KV = L 

€F < AKSIKA ? ILi 12) ® EQ^ INITALIKV)) GO TO 30 

29 CONTINUE 

30 KN ( KV ) =0 

t F t KV .N€. Kl) GO TO 34 
DO 32 KK*l t N 
Kt-KMIKKt 

■IF ( K L ® NE ® 01 GO TO 33 

32 CONTINUE 

IFlOyTUU It 12) *NE« CM) GO TO 360 
0UT<J»IUI2)±RP 
JRP=sJ 
GO TG 36 

33 K l-KL 
K A - KL 
GO TO 35 

34 IFICUTI J*U 9 i2).EQ. LP) GO TO 335 
J*= J + l 

OUTI Jt-I H L2H AN ST KA| 1 1 f 1 2 ) 

532 DO 332 KK=1,N 
KL “KN { KK i 

1FIKL .N€. 0) GO TO 333 

332 CONTINUE 

333 IF (KL ® NE o Kl) GO TO 334 
KN (KL)=0 
GO TO 532 

334 4FIK1 .EC. 01 GO TO 435 
J = J41 

OUT Id *1 Mr 12 ) ^CM 
KA = KV 
GO 10 35 

435 

QUTI J, I If I2)^RP 
URP^d 
GO TO 36 

335 JfcJTl 

0UT(J»I UI2)MNI TALIK1) 
id - J T 1 

mi itif l IV I 2 ) -CM 

OUTI Jf I u I2)±AN$(KA f f l t 12) 



o o 


<F(A-NS(K'Vt 11(121 • EQk INITAL (KA) ) GO TO 336 
J«J*1 

OUTI Jt-1 l f I2)±CM 
&A~K V 
GO TC 35 
336 KN { K A )~0 
3-J+l 

0UTf J„IU I2MRP 
632 DO A 3 2 KK ~ l ? N 
KL-KNIKK1 

I F { Kl_ • N€ * 01 GO TO 433 

432 CONTINUE 

433 I F { K L .N£. K 2 ) GO TO 434 
KN ( KL 0 

GO TO 632 

434 *)RP-0 

I F { KL .EC® 0) GO TO 36 

0011 JUli I2)^LP 

K1»KL 

&A = KL 

35 CONTINUE 

360 IFIGOT( JU It 12) e NE« LP) GO TO 36 
OUT OJ f I If I2)^BLK 

36 J1=J 
J^JRP 

€ F ( J ® NE ® € ) GO TO 37 
OUT (J 14 lit I2)=BLK 

37 NOO Till , 12 )=d 

IF (U*LE.NAX! GO TO 38 
f^AX~ J 

3E CONTINUE 

BLOCK C 7 ) DECODES AND WRITES OUTPUT 

DO 42 L-lf576 
42 L INE I L ) =BLK 

€KY = C 

DO 45 L 2 fc 1 1 2 

NCHK^C 

NPER^O 

DO 45 L 1*1*400 

1FINPER ® GE ® NOPAIL23) GO TO 45 
€ F < IfcPUT<Ll,L2) ®EQ® PER) GO TO 445 
NCHK&NCHK^l 
MCFKlsNCFK 
GO TO 45 
44 5 NL = L 1 — NC F K— 1 
NPER^NPER^l 
DO 43 I -l t NCHK 
NL=NL+1 

L INPUT ( I i =5 INPUT f NL f L2 ) 

43 CONTINUE 
LT AG S 1 

CALL RE AG C L I N PUT ? L I NE f L T AG ) 

NCHK-C 
45 CONTINUE 

00 101 I * 1 1 NO PI 
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1-01 4TCPCI)=tlNEi I) 
iNC s C 

NCR HI 

NAPA1HN0P 1 + NGP2 
00 i 02 1 ~N A * NAM AX 
IND^IND+i 

10 2 IS ID I INDI-LINEU ) 

DO 46 L s lf NAMAX 
46 LINE ( L ) -8LK 
IK Y = C 
LTAG=C 
NCHK* 15 

DO 52 KO fe 1 , N0P2 
DO 52 JO^lfNGPl 
DO 50 10*1 f 15 

50 L. INPUT ( IG ) =GUT! I 0 *JOi KO ) 

CALL RE AC ! L INPUT ? L I N£ , L T AG ) 

52 CONTINUE 

WRITE 16^68) TITLE 

KSN = C 

LSW=€ 

DO 115 1L = H6 
IL1* IL-l 

44= ( KSW.EC.l) GO TO LI 5 
MAX 3 20^ I L 

IFINCP1 ® GT o 119) GO TO 111 
IF ( NOP 1 wGT * MAX) GO TO 104 
MAX 3 NOP 1 
KSW*1 

€ F ( NCP 1 sLE, 20) MAX* N0P1 

104 1N~0 

ITEM* 20*HL1*1 
DO 1C5 I J=? I TEM ? M AX 
IN* IN+ 1 

105 OUTPUT! I N ) = I TOP ( I J) 

INI 3 IN + 1 

DO 1C55 IP=INl,21 

ig=ip+i 

TD0T4 LG ) fe BLK 

1055 OUTPUT! IP)= BLK 

1056 WRITE (6,11041 < OUTPUT U) , I = 1 ,20 1 
1104 FORMAT! IFG,1QX,20A6) 

WRITE (6 ,1 144 1 { I DOT { T ), 1*1,21) 

1144 FORMAT! 1H , 3X,-21 A6) 

DO 1145 1 R-l 1 21 
114 5 I DOT I IR)*IDIT 
LB*C 

106 00 1C 8 LN 3 1 , N0P2 
LG*1 

LANN*LB4LN 

I F ! L ANN *GT® 119 ) GO TO 109 
OUTPUT ( H * I SI D(L ANN ) 

DO 107 I J* 1 , MAX 
NMAX*20*ILH! J 
f F ! NW AX kGl« NOP 1 ) GO TO 107 
LG=LG* 1 

IfUi *LT* 21 GO TO 1076 

MAR* 20*! 1 L+LN-2 ) ♦ C NOP H 20 ) * ( LN™ 1 ) + IJ 

SO TO 1077 



l§76 WAR=?NGP1M It +LN-2 > + f J 
107 7 OUTPUT ( LG ) =L I NE C MAR ) 

107 CONTINUE 
LG = LG+'l 

IFCLG «GT« 21) GO TO 1089 
00 1 C 88 TQ - LGt 21 

1068 OUTPUT ( IC)=BLK 

1069 CONTINUE 

108 WRIT€<6, 11051 ( OUTPUT (I ) U = 1 f 2 1 1 
1105 FORMAT! 1H * 4 X f A6 * 20A6 ) 

109 CONTINUE 

T F C L S W *FQ. 1) GO TO 115 
I F ( N G P 2 uL£. 56) GO TO 115 
IFIN0P2 *GT« 112) GO TO 110 
L S K = 1 
LB = 5 6 


GO TO 106 

110 IFCNCP2 . 01 . 120) GO TO 111 
L'Sh-1 

LB^ 112 
GO TO 106 

111 WR I TE < 6 1 1 1 06 ! N0P1 ? N0P2 

1106 FORMAT ( 1FK,39HT00 MANY ROW OR COLUMN OPERATIONS NOP 1= I4§ 5HN0P2- 14 I 
115 CONTINUE 
GO TO 1 

55 WRITE (6*711 ( IOENT ( I ) * IM* ID) 

WR I TE ( 6 t ? 9 ) TITLE 

80 TO 1 

56 WRITE ( 6*72) 

GO TO 1 

57 WRITE (6*73) 

GO TO 1 

58 WRITE (6% 74) 

GO TO 1 

59 WRITE (6*75) 

GO TO 1 

60 WRITE ( 6*76) 

GO T € 1 

61 WRITE (6*77) 

GO TO 65 

62 WRITE (6*72) 

GO TC 65 

63 WRITE (6*73) 

GO TO 65 

65 RETURN 

66 FORMAT < 80 A 1 1 


68 FORMAT! 1HK,5GX,20A1 ) 

71 FORMAT! 1H J * 44HI L LEGAL ELEMENT IN GROUP. IDENTITY GROUP IS , 11A6) 

72 FORMAT ( IF J * 29HBL ANK IS NOT A VALID ELEMENT ® ) 

73 FORMAT ( IF J ? 5 1HILLEGAL USE OF PARENS OR TITLE + IDENT CAROS MIXED) 

74 FORMAT ( IP J ? 14HI L LEGAL GROUP. I 

75 FORMAT! IF J ? 34HTHE PAIRS ARRAY HAS BEEN EXCEEDED. ) 

76 FORMAT! !HJ f 34HTOO MANY OPERATIONS. LIMIT IS 24.) 

77 FORMAT! lFJ f 39HM0RE THAN 6 ELEMENTS IN IDENTITY GROUP.) 

79 FORMAT! 1FJ?27HCHECK INPUT CARDS IN GROUP f 80A1 ) 

80 FORMAT I 1F1 f 80A1 1 

81 FORMAT! 1HJ # 80A1! 

9 C FORMAT I IT J f 80Al) 

END 
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$IBFTC READ 


C THIS SUBROUTINE CAN ACCOMODATE BOTH S3 AND S4 

C 

SUBROUTINE R E AD I I NPUT ♦ L I NE » L T AG ) 

DIMENSION INPUT { 15) ,LINE<576), I GROUP ( 6 ) 

■INTEGER SRl,SR2, SR3 , SR4 » SR l S ,SR2S , SR3S , SR4S , R 1 , R2 » R 3 , R 1 S , R2S , R3S , 
.R1C »R 2C , R 3C 

INTEGER ONE, TWO, THREE, LP, RP 
DATA IDNT , LP i RP/ 1HE,1H( , 1H)/ 

DAT A I IGRCUP ( I ) , I = 1 , 3 I / 1H1 , 1H2 , 1H3/ 

COMMON/ IKYMAX/IKY,MAX»NCHK 
COMMON/ NO MOPS/JO, KO,N OUT I 24, 24) 

CO M, MON/LI TRAL/I I 1,112,13,14,15, 1 6 , SR 1 , SR2 , SR3 , SR4, SR1S , SR2 S , SR3S , 
«SR4S ^ R1 ,R2,R3,R1S,R2S,R3S,R1C,R2C,R3C 
0NE=f GROUP S II 
.TW0=4 GROUP (21 
THREE=I GROUP! 3) 

IFSLTAG » EG . 1) GO TO 1 
IFINCUTI JO, KG) . NE. 0) GO TO 1 
>IKY = IKY+1 
LINE! IKYI= IDNT 
GO TC 19C0 
1 DO 18 CO 1=1 , NCHK 

•IF (INPUT! I) .EQ. LP) GO TO 1800 
IF ( INPUT I I ) .EQ. ONE) GO TO 1100 

I F ( I NPUT ( I ) .EQ. TWO) GO TO 1010 

3 F ( INPUT 1 1 ) .NE. THREE! GO TO 1800 
4KY=f KY+1 
LINES IK Y 5 = 1 1 2 
GO TC 19C0 

1010 4 F ( INPUTS 1+2) .EQ. THREE) GO TO 1020 

I F 1 1 NPUT ( I +3 1 .EQ. RP) GO TO 1011 

IKY=IKY+1 

LINE S IK Y) 3 SRI 
60 TO 1 9€0 

1011 4KYMKY + 1 
LINE1 IKYHI5 
GO TC 19C0 

1020 I F ( INPUT I I +3) .EQ. RP) GO TO 1021 
IK Y = IKY+,1 

LINES IKY)=SR1S 
GO TC 1 9C0 

1021 IKY=f KY+1 
LINES IK Y ) =5 1 4 
GO TC 1 9C0 


1100 1 F ( INPUT I 1+21 

.EQ. TWO) 

GO 

TO 

1200 

IF ( INPUT! 1+2) 

. EQ® THREE ) 

GO 

TO 1130 

IF! INPUT! 1+3) 

® NE® RP) 

GO 

TO 

1110 

IF ( INPUT 11+41 

a EQ® LP) 

GO 

TO 

1101 


IKY=IKY+1 
LINES IK Y)=I3 
GO TG 1 9C0 
1101 1 K Y = I KY +,1 

LINES IKY) = R3S 
GO TC 19C0 

1110 3 F ( INPUT ( I +4 ) .EQ. TWO) GO TO 1120 
IF ( INPUTS I +5) .EQ. RP ) GO TO 1111 
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mi 
112 0 

1121 

1130 

1131 

1140 

1141 
115 0 

1151 

1200 

1201 

1210 

1211 

1120 

1221 


iKY~f KY*1 

LINE! LKYMR1C 

GO TC 1 9C0 

IKY- IKY*1 

LINE! IKY} = SR4S 

€0 10 1 9€0 

•I F ( INPUT! 145! .EG* 

4K Y= IKY 4.1 

LINE I IK Y ! = R2 

GO TO 19C0 

4 K Y = I K Y 4 1 

LINE! IKY}=$R3 

GO TC 1 9€ 0 

IF ! INPUT I I 43) „NE« 

IF! INPUTII441 oEQ. 

IKY=IKY4l 

LINE! I K Y ) =j I 6 

GO TC 1 9C0 

!KY=f KY41 

LINE! IKYI=?R1S 

GO TG 1 9€0 

IF( INPUT! 1441 *EQ* 

IF ( I NPUT 1 1 45 ! ,EQ, 

1 KY = I KY 4 1 

LINE! LK Y )=R3 

GO TC 1 9C0 

4 K Y = 1KY41 

LINE! IK Y 1 =?SR4 

GO TG 1 9C0 

1 F { INPUT I 145} «EQ« 

I K Y= I KY 4.1 
LINE! IK Y > = R2C 
GO TG 1900 

lKY=rIKY4l 

L INE ! IKY I =? SR2S 

GO TC 19C0 

IF ! INPUT! 143} *NE* 

IF (INPUT! 1441 «EQ® 

I KY = I KY +1 

LINE ( IKYI =11 

LINE! IKY}= III 

60 TC 19C0 

4KY = i KY+.-l 

LINE I LKY) = R2S 

GO TO 1900 

IF (INPUT! 1441 *EQ. 

IF! INPUT I 1451 .EQ. 

IK Y= IKY 41 

LINE ! I K Y1 = R 3C 

GO TG 1900 

IK Y=i KY +1 

LINE! I K Y } = SR 3 S 

60 TG 1900 

IF (INPUT! I 451 .EQ. 

IKY=4KY4l 

LINE I I K Y ! = R 1 
GO TG 1900 
4 K Y = l KY 4.1 
LINE i IK Y J - SR2 


RP) GO TO 1121 


RP) GO TO 1140 
LP) GO TO 1131 


TWO! GO TO 1150 
RP) GO TO 1141 


RP) GO TO 1151 


RP) GO TO 1210 
LP) GO TO 1201 


THREE) GO TO 1220 
RP) GO TO 1211 


RP) GO TO 1221 
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GO TO 1 9€0 
1800 CONTINUE 
1900 RETURN 
END 


APPENDIX C 



2 


















IKY = IKY + 1 
LINE(IKY) = R2C = (1, 3, 2, 4) 
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